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Abstract 

To investigate the characteristics of atmospheric gasoil as a 
suitable feed stock of steam cracker processes, several 
experiments were designed and carried out in a pilot plant. 
The operating conditions of the designed experiments were 
coil outlet temperature (COT), feed flow rate and steam ratio 
(STR). A group of these experiments were conducted to 
determine the coke formation rate. Central composite design 
methodology was utilized to set the experiments and 
analysis the results. Some statistical tests were used for 
confirmation of the accuracy, consistency and reproducibility 
of experimental results. Moreover, statistical model was 
utilized for investigation of the coke formation rate. To 
determine the yield distribution of main products a reaction 
network consists of 21 reactions together with the related 
kinetic reactions rates were developed. Based on the 
developed models and utilizing suitable optimization 
algorithm, the best operating conditions were determined. 
The results declared that at the 883.9°C as coil outlet 
temperature,0.62 as steam ratio and residence time equal to 
0.3 sec the optima levels were obtained in which the yield of 
ethylene, propylene and butadienes were respectively equal 
to 28.73%,12.33% and 4.02%. 
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Introduction 

Olefins are widely used as basic species in production 
of polymeric material in petrochemical complexes. 
Ethylene as one of the main light olefins is widely 
produced in the world. It's globally production was 
107 million tones in 2005 and 109 million tones in 
2006.By 2010 ethylene was produced by at least 117 
companies in 55 countries. To meet the ever increasing 
demand for ethylene, sharp increases in production 
facilities are added globally, particularly in the Persian 
Gulf countries and in China [l].Over 80% of ethylene is 
used to produce ethylene oxide, ethylene dichloride 
and different polyethylenes. Between products, low 
density and linear poly ethylene are widely demanded 



products [2]. Propylene as the other light olefins has 
been mass produced since sixty years ago. It is used to 
produce propylene oxide, polypropylenes, propanol 
and the other related chemicals [3]. Steam cracking 
process is the main commercial unit for mass 
production of these olefin species in petrochemical 
complexes. In this process, gas or liquid hydrocarbons 
are heated to temperature in the range of 770 °C to 
870°C .In this range of temperature, the radical 
reactions are accelerated. These kinds of reactions 
convert heavy hydrocarbons to lighter ones and the 
saturated species are converted to unsaturated. The 
complex mixture is yielded from the proposed reaction 
network which contains propylene and ethylene. The 
produced mixture is sent to some separation 
equipments to separate and purify the valuable 
produced species. 

Coke as one of side product of the process is produced 
in the secondary steam cracking reactions and settled 
down in the reactor and TLE (Transfer Line 
Exchanger). This reduces the heat transfer rate, 
increases the pressure drop and the possibility of 
formation of hot spots. The amount of formed coke 
depends on the feed stock characteristics, operating 
conditions and the alloy used in the reactor pipe. The 
amount of contaminated coke on the inner wall of the 
pipe is increased by passing the time. To lower the 
effects of coke on the efficiency of the process the 
reactor temperature are increased. This also increases 
the coke formation rate and at last it causes to shut 
down the reactor furnace to clean up [4]. 

Coke formation is a complex phenomenon which 
follows from four different mechanisms: 

1. Catalytic coke formation 

2. Pyrolytic coke formation 

3. Condensation mechanism of coke formation 

4. Transformation of coke to down streams 
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In the first step, the heterogeneous reactions are 
carried out on the metal surface as heterogeneous 
catalyst of reactor pipes and consequently the 
filamentous coke is formed [5]. This has the highest 
rate in comparison with the other coke formation 
mechanisms while the operating temperature is low 
and in the range of 500°C to 600°C and the process is at 
the starting levels. In the pyrolytic coke formation 
mechanism which is also known as radical coke 
mechanism, the coke precursors joint to the coke 
surface and grow the coke [6]. In this mechanism 
several species as acetylenes, olefins, aromatics and 
paraffins incorporate in different reactions to form this 
kind of coke. The highest reaction rates belong to the 
reactions in which the acetylene is involved and the 
lowest is for the reactions with paraffins as reactants. 
The produced coke has amorphous structure [7]. The 
other steps of coke formation as condensation 
mechanism and coke transformation are highlighted at 
the end of reactor while the operating temperature is 
decreased. 

Due to vitality of coke formation in thermal cracking 
processes, several researchers focused on this subject. 
The effect of operating conditions as residence time, 
temperature and partial pressure are investigated by 
Kunzru and Kumar [8]. Sundram and Froment [9] 
conducted other research to develop a rigorous model 
to predict the coke formation. In the proposed model 
the propylene used as the major coke precursor. The 
effects of alloys and materials used in the reactor pipes 
and the coke precursors in gas thermal cracking are 
studied by Zou et al [10]. Marin & Wauters [11], 
developed a more precise model to investigate coke 
formation rate. The proposed reaction network had 
14000 reactions and 2400 species for ethane thermal 
cracking. 

To inhibit or to lessen the coke formation, several 
researchers focused on the coke inhibitors [12-15]. The 
effect of DMDS as one of the coke inhibitors are 
studied by Dhuyvetter et al [12] .Moreover, the effect of 
IHbPtCk as coke inhibitor is studied by Chan et al 
[13]. The power of Thiophene and benzothiozol to 
inhibit the formation of coke in thermal cracking of 
atmospheric gasoil is clarified in the Shubo et al 
research [14]. Also, the effect of a mixture of 
components contained Zn, Si and Sulphurous 
compounds is investigated by Brown et al [15]. 

In this research, by carrying out several experiments 
designed and analyzed based on central composite 
design methodology the effects of different operating 



variables on coke formation are estimated. Also, a 
statistical and a rigorous kinetic model are developed 
to predict the coke formation rate and the yield of 
main primary and secondary products. Then the 
optimization is carried out to determine the optimal 
point and the best reactor temperature profile. 

Experimental Section 

Feed Stock Specification 

The feed stock of steam cracking pilot plant is an 
atmospheric gasoil with a boiling range of 218 Q C to 
387 Q C, with a density of 0.845 gr/cm3. The total sulfur 
of the feed stock is approximated around 0.75 wt%. 
The main specifications of this oil cut are listed in 
Tablel. 

TABLE 1 THE MAIN PROPERTIES OF THE ATMOSPHERIC 
GASOIL 



Specification 


Gas-oil 


Nitrogen, Total wt% 


<0.5 


Hydrogen, Total wt% 


13.8 


Carbon, Total wt% 


85.6 


Aromatic Content, vol% 


14 


Olefin content, vol% 


Trace<0.3 


Saturate Content, vol% 


86 


Distillation (D86): 


°C 


IBP at 760 mm Hg 


218 


5%Vol. Recovery 


257 


10%Vol. Recovery 


266 


FBP at 760 mm Hg 


387 



Pilot Plant Setup: 

The selected experiments were performed in a steam 
cracking pilot setup that was designed and assembled 
for steam cracking of the hydrocarbon feed stocks in 
the range from ethane to vacuum gasoils. The 
hydrocarbon feed stock and the diluted water are first 
preheated to the hydrocarbon cracking 
temperature(approximately 550 Q C) then mixed 
together and fed to the reactor in which the main 
cracking reactions are take place and the reactant 
mixture are heated up to 860-C. The reactor effluent 
would be immediately quenched by cooling water in a 
double pipe heat exchanger. In order to separate the 
condensate from the gaseous product stream, the 
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exchanger outlet stream is sent to a flash drum. Then, 
the gas phase enters a series of condensers, passing 
through a filter in the final stage. 

Experiments 

Central composite design methodology was applied to 
design the experiments. This is used with three design 
factor, namely, the hydrocarbon feed flow rate (XI), the 
steam ratio(X2) and coil outlet temperature (X3). The 
coded level and the natural values of the mentioned 
variable are clarified in Table 2. 



The number of trials was based on the number of 
design factors and was equal to 19 experiments (15 
combinations with four replications). The results of 
experiments and the design matrix are shown in table3. 

As introduced in Table3, the maximum coke formation 
rate is obtained while the coil outlet temperature, 
steam ratio and feed flow rate are respectively equal to 
1,-1 and -1 in coded level, while minimum coke 
formation rate is obtained for coil outlet temperature, 
steam ratio and feed flow rate equal to -1,1 and 1 in 
coded levels. 



TABLE 2 CODED AND NATURAL LEVELS OF THE DESIGN FACTORS. 



Design Factors 


-1.6818 


-1 





1 


1.6818 


Feed Flow rate 


0.98 


2 


3.5 


5 


6.00 


Steam ratio 


0.4636 


0.6 


0.8 


1 


1.13636 


Coil Outlet temperature 


716 


750 


800 


850 


884 



TABLE 3 DESIGN MATRIX AND RESULTS OF THE CENTRAL COMPOSITE DESIGN 



No 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


15 


16 


17 


18 


19 


Feed Flow 
rate 


-1 


-1 


-1 


-1 


1 


1 


1 


1 


-1.68 


1.68 























-1 


-1 


Steam 
ratio 


-1 


-1 


1 


1 


-1 


-1 


1 


1 








-1.68 


1.68 

















1 


1 


Coil Outlet 
temperature 


-1 


1 


-1 


1 


-1 


1 


-1 


1 














-1.68 


1.68 











1 


1 


Rate of 








































Coke 
(gr/cm 2 .sec) 
10 7 


1.67 


34.1 


0.66 


20.9 


0.16 


18.2 


0.03 


3.9 


12.3 


0.51 


11.53 


0.65 


0.06 


32.4 


1.6 


2.35 


2.4 


24.6 


25.6 



Modeling Section 

Two different models as statistical and rigorous kinetic 
model are used to investigate the coke formation rate 
and product yield distribution. The following full 
quadratic model was used to develop the statistical 
model. It is obtained by a multiple regression 
technique for three factors. 

Yi=a+[3iXi+[3 2 X2+ (3 3 X 3 + [3nXi 2 + [3 22 X2 2 + [33 3 X 3 2 + |3i 2 Xi X 2 + 
|3i 3 Xi X 3 + [3 23 X 2 X 3 (1) 

In table 4, the significance of different factors and their 
interactions, the related coefficient, the ANOVA 
results for lack-of-fit and the result of R-squared tests 
have been shown. 

This table introduces the significant and insignificant 
parameters by declaring the p- value. The vital value 
for "p" is 0.05. The values less than this border point 
introduce the significant parameters. The insignificant 



parameters are discarded with a confidence level of 
95% while their p-value are higher than 0.05.The p- 
value of lack of fit is used to determine whether or not 
the constructed model was appropriate to describe the 
observed data. 

To investigate the yield distribution of main products 
a molecular reaction network is introduced. Based on 
the introduced reaction network, reaction mixture is 
divided into two different sections as primary 
pyrolysis and dissociation sections. The primary 
reactions are the only chemical phenomena happened 
in the first section. However, in the second section the 
decomposition and polymerization reactions occurred 
and the primary and secondary reaction product 
accompanied as the reactants. Pyrolysis-oil which 
contains several aromatic species is one of the second 
section reactions products. 
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TABLE 4TEST OF SIGNIFICANCE OF FACTORS AND 
INTERACTIONS FOR STATISTICAL MODEL OF COKE 
FORMATION, R-SQUARED AND LACK-OF-FIT 



Factor or 
Interaction 


Rate of Coke formation * 10 7 
(gr/cm2.sec) 


Coefficient 


T 


p-value 


Constant 


968.57 


7.231 





Xi 


38.013 


5.41 





Xi 


193.05 


3.796 


0.004 


x 3 


-2.92 


-8.838 





Xi*Xi 


0.47689 


2.093 


0.066 


x 2 *x 2 


16.388 


2.438 


0.037 


x 3 *x 3 


0.002193 


10.442 





Xl*X2 


-0.6494 


-0.294 


0.775 


Xl*X3 


-0.05453 


-6.517 





x 2 *x 3 


-0.28708 


-4.795 


0.001 


R-Squared 


99.10% 


Value of lack of 
fit 


0.073 



To demonstrate the primary reactions, the following 
pseudo reaction is introduced [16]: 



Feedstock - 



±->aH 2 + bCH 4 + cC 2 H 4 + 



dqH 6+ eC 3 H 8+ fC 3 H 6 + 

g QH 8 + h QH 6 +iC 2 H 2 + jC 4 H 10 



(2) 



Second part of reactions network contain several 
reactions to decompose or combine the ten products 
yielded from reaction demonstrated in 
equation(2).Table(5) clarifies reactions happens in this 
section. In this table, EB and B refer to Ethyl benzene 
and benzene respectively. To consider, other aromatic 
species a pseudo species "C6.5 H7"is defined. Also "L" 
refers to pseudo component with 7.5 carbons. 



TABLE 5. THE REACTIONS OF SECONDARY SECTIONS IN REACTION NETWORK 



Reaction 

1MO 


Reaction 




■(") 

\mol J 


1 
1 


C H C H -1-/7 

2 6 2 4 2 




Z/o 


2 


ZL -2 n 6 *^3 n 8 + LU 4 


1.87xl0 n 


275 


3 


C 2 H 4 + C 2 H 6 >C 3 H 6 + CH 4 


2.57xlO n 


253 


4 


C^H g <-> C^H ^ + H 2 


3xl0 10 


214.7 


5 


^3^8 ^C 2 H 4 + CH 4 


1.03xl0 12 


211.8 


6 


C 3 H 8 + C 2 H 4 >C 2 H 6 + C 3 H 6 


1.267xl0 13 


247.2 


7 


2C 3 H 6 >3C 2 H 4 


1.636xl0 12 


268.7 


8 


2C 3 H 6 >0.3C 65 H 7 + 0.14L + 0.3CH 4 


3.6xl0 10 


208 


9 


C 4 H W ^2C 2 H 4 +H 2 


3.5xl0 16 


295.9 


10 


C 4 H 10 >C 2 H 4 + C 2 H 6 


1.035xl0 15 


256.6 


11 


C 4 H W < >C 4 H 8 + H 2 


1.53xl0 14 


261 


12 


C 4 H 8 >0.41C 65 H 7 + 0.19L 


1.432xl0 5 


212.4 


13 


C 4 H 8 >C 4 H 6 +H 2 


1.182X10 11 


209.3 


14 


C 2 H 4 +C 4 H 6 >B + 2H 2 


8.3xl0 2 


103.2 


15 


C 4 H 8+ C 3 H 6 ^T + 2H 2 


1.38xl0 4 


140.4 


16 


C 4 H 6 + C 4 H 8 > EB + 2H 2 


3.613xl0 12 


246.1 


17 


2C 4 H 6 >ST + 2H 2 


1.32xl0 12 


263.4 


18 


C 2 H 4 >0.15C 2 H 6 +0.233C 3 H 6 +0.1C 4 H 8 


3.148xl0 15 


303.8 


19 


C 3 H 6 + C 2 H 6 >C 4 H 8 + CH 4 


9.673xl0 13 


291.1 


20 


C 3 H 6 >C 2 H 2 + CH 4 


9 


59.38 
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To determine rate constants as frequency factor and 
activation energy, least square algorithm was utilized. 
The objective function clarifies in the equation (3) is 
used in the algorithm: 

/ = W l X (Cjmodel ~ C jexp ) 2 + W 2 £ (T jmodd - r. exp ) 2 (3) 

i 

Cj model and Cjexp are respectively the reactor outlet 
molar concentration calculated from model and 
obtained from experiments. Moreover, Tjmodel and 
Tjexp are the temperature of different reactor sections 
fixed or determined in model and experiments. The 
procedure of investigation of rate constants is 
demonstrated in figure (1). 



As demonstrated in figure 1, the first step of the 
procedure is generation the initial estimations. These 
are the initial estimates of rate constants and activation 
energies of the proposed reaction network .In the 
second stage, the reactor model is run with the 
introduced constants to calculate the outlet 
concentration of different species. In the next step, 
utilizing the calculated and experimental information 
the objective functions are investigated. If the objective 
functions values meet the minimum value, they are 
reported as the optimum value. Otherwise, the next 
generation of rate constants are generated and fed to 
the other steps of the determination procedure. 




Initial Estimate generation 



Modeling Calculation 



Optimization 



Objective function 



No 




FIG. 1 THE OPTIMIZATION PROCEDURE FOR DETERMINATION OF KINETIC RATE CONSTANTS 
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Result and Discussion 

Based on developed statistical model, the effects of 
variation of operating variables on coke formation rate 
are determined and demonstrated in figure(2) to 
figure(4).In figure(2), the effect of variation of feed 
flow rate and coil outlet temperature(COT) is clarified. 
Increasing COT, less than 760(°C) for stable feed flow 
rate decreases the coke formation rate. On the other 
hand, an increment in feed flow rate for stable COT 
decreases the coke formation rate. 




CQTiC) 




FIG (2) - COKE FORMATION RATE VERSUS FEED FLOW RATE 
(GR/MIN) AND COIL OUTLET TEMPERATURE (°C) AT 
CONSTANT STEAM RATIO (GR/GR) 

As figure (3) declares, increasing steam ratio with 
constant operating conditions decreases the coke 
formation rate. Also, equal increment in steam ratio 
and COT raises the rate of coke formation with the 
style like as the position in which only the temperature 
is increased. 



1.1 
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/ ///// 


m 


07 
0.6 
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FIG (3) - COKE FORMATION RATE VERSUS COT (°C) AND 
STEAM RATIO AT CONSTANT FEED FLOW RATE (GR/MIN) 

Moreover interaction effect of feed flow rate and steam 
ratio on coke formation rate can be obtained in figure 
(4). Indeed, equal increasing of steam ratio and feed 
flow rate decreases the rate of coke formation. 




FIG (4) - COKE FORMATION RATE VERSUS STEAM RATIO 
(GR/GR) AND FEED FLOW RATE (GR/MIN) AT CONSTANT COIL 
OUTLET TEMPERATURE (°C) 

The maximum observable coke formation rate which is 
equal to 6.62*10 6 gr/cm 2 .sec is observed at the 
maximum coil outlet temperature, minimum steam 
ratio and minimum feed flow rate which are 
respectively equal to 884( Q C), 0.464 (gr/gr) and 0.977 
gr/min. 
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FIG (5) - YIELD OF DIFFERENT PRODUCT VERSUS YIELD OF 
METHANE 

As mentioned above, a molecular reaction network 
was developed for modeling of steam cracking 
reactors. Figure (5), demonstrates the effects of the 
variation of methane yield on the yield of different 
products. Increasing the yield of methane always 
decreases the yield of C5+ and increases the ethylene 
yield. However, the propylene yield has limited 
variation by increasing the methane yield. Applying 
the developed statistical and kinetic models together 
with SQP (Sequential Quadratic Programming) 
algorithm as a powerful optimization tool, 
optimization of the pilot plant is conducted. 
Maximization of net profit of the process is usually 
suitable objective function in optimization procedure 
of a chemical process which can be defined as follow: 



/ = income - cos t 



(4) 



In the above equation, the "income" refers to the price 
of products. Moreover, "Cost" refers to the production 
cost, in which the cost of coke regeneration, energy 
consumption during steam cracking and the cost of 
consumed water are included in production cost. 

Hence, the "income" item mentioned in equation (2) 
should be introduced as follow: 

income = 

6 * 10 4 * mfeed * (1 . 15 * Ethylene Yield + 
1.16* Propylene Yield + 0.831* Butadiene Yield + 
0.35 * ButenesYield + 0.266 * Fueloil Yield + 0.575 * Lightgas Yield) 

(5) 

The cost of process refers to the cost of consumed 
water, cost of feed, cost of heat and consumed 
energy .Cost of water and feed flow rate are 
encountered in the following relations: 



o o 
cos t = 2.64 * 10" 4 * m feed * steam _ ratio + 0.03 156 * m feed 

(6) 

TABLE 6 PRODUCT VALUED AND OPERATING COSTS [17, 18] 



Item 


Value 


i ower^/Kvvrij 


u.uo 


reea^/Kgj 


U.DZD 


/vcetyiene^/Kgj 


\J.yD 




1.15 


Propylene($/kg) 


1.16 


Butadienes($/kg) 


0.352 


Deionized water($/kg) 


0.0044 


Fuel oil($/kg) 


0.266 


Light gases($/kg) 


0.574 



The cost of consumed energy is divided into two parts 
as the consumed energy for olefin production and 
consumed energy for decoking process. The price of 
consumed energy for production of olefins is 
calculated by the following equation: 



0.216 
1 



* UA *( T skin- T ) 



(7) 



In this equation, r), TSkin and T respectively refer to 
the furnace efficiency, skin temperature of reactor pipe 
and the temperature of reactive mixture. The cost of 
decoking is calculated by equation (8); 



P = P + P 

1 2 1 2C ^ 1 2F 



(8) 



P2C is the cost of compressor consumed energy; P2F is 
the cost of energy consumption in furnace during the 
decoking procedure. The consumed energy is 
calculated by the following equation. 



p —(^ r )^f t *n *p * ( Met \ v * 

r 2c — w.^iu \J Mr r outlet v ) 



k P — 

[( ouotlet ) k " 



k-1 



The P2F variable is defined as follow: 



p =p * 2 

r 2F r l /l 



(9) 



(10) 



A is the ratio of decoking over olefin production time 
periods. The defined optimization problem has some 
constrained. At first, due to the metallurgical 
limitation, the reactor temperature should not be 
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exceeded than 890-C. As the second constrained, the 
yield of each product have to be in the range of 0% to 
100%. And due to the mass continuity rule, the reactor 
mass inlet flow rate should be equal to the reactor 
mass outlet flow rate. The result of optimization 
procedure is calculated and clarified in table7. 

TABLE 7 THE RESULT OF OPTIMIZATION 



Operating condition 


Value 


Coil outlet Temperature( Q C) 


883 


Feed flow rate (gr/min) 


6.02 


Steam ratio(gr/gr) 


0.62 


Residence time(sec) 


0.3 


Y Acetylene(%) 


1.23 


Yethylene(%) 


28.23 


Ypropylene(%) 


12.33 


YButadiene(%) 


4.02 


YButenes(%) 


0.98 


Ylg(%) 


9.70 


Yfueloil(%) 


28.16 


Rate of Coke formation 
(gr/cm2.sec)*10 7 


24.32 



Determination of optimum temperature profile in 
reactor is one of great challenges in running of steam 
cracking process. Three different temperature profiles 
for three different optimization cases as maximization 
of ethylene production, maximization of propylene 
production and maximization of profit are determined 
and demonstrated in figure(6). 



Ethylene maximization 



Profit maximization 




FIG. 6THE OPTIMUM PROFILE FOR MAXIMIZATION OF PROFIT, 
PROPYLENE AND ETHYLENE AT RESIDENCE TIME AND 
STEAM RATIO RESPECTIVELY EQUAL TO 0.362 SEC AND 0.88 

As demonstrated in the aforementioned figure, 
Temperature profiles for two optimization cases which 
are solved to maximize the yield of propylene and 
ethylene have many similarities. The reactor inlet and 



outlet temperature are same in both profiles. But, the 
average reactive mixture temperature for ethylene 
maximization case is greater than for the propylene 
maximization challenge. Infact, the rate of temperature 
increment for the first case is greater than for 
temperature profile of the second case. The maximum 
point of temperature profile of the second case is 
greater than the first one. The third one has many 
differences from the others. The style of these profiles 
can be interpreted by the nature of production of the 
two important olefins as ethylene and propylene. 
These two olefins are the main primary products 
which contribute in formation of the secondary 
products. But propylene has greater role and 
consequently is more consumed in secondary products. 
Increment the temperature along the reactor with high 
rate, raises the possibilities of improvement of 
secondary reactions that will cause to consumed more 
propylene. So, reduction of the rate of temperature 
raising at the inlet reactor regions and maximizing the 
temperature before sharp declining maximized the 
yield of propylene. But, since ethylene is produced in 
primary and some secondary reactions, determination 
of the operating condition to make little improvement 
in the secondary reactions may maximize the ethylene 
production. Moreover, the temperature profile for 
profit maximization shows several differences from 
the others. The inlet, average value and the outlet 
temperatures of this profile are less than the others. 
Infact, the demonstrated profiles declares that the 
amount of energy consumption has highlighted effect 
in the profit of the process. 

Conclusions 

To evaluate the atmospheric gasoil as a suitable feed 
stock of steam cracking process, several tests were 
carried out in a pilot plant. Based on the experimental 
results, a rigorous kinetic and a statistical model were 
developed. The results of statistical model proposed 
that increasing the Temperature could increases the 
rate of coke formation. While increasing the feed flow 
rate and steam ratio has inverse effect and could 
decrease the coke formation rate. 

Based on rigorous kinetic model and the experimental 
results, the yield of main products as ethylene, 
propylene and butadienes versus the yield of methane 
as an index of improvement of reaction network were 
studied. While the yield of ethylene shows increasing 
trends, the style of Cs+ yield has decreasing trend and 
propylene has limited variation. The optimization was 
carried out to determine the best operating conditions 
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for profit maximization. To develop the related 
objective function, the price of product is considered as 
income and the summation of the price of feed, cost of 
consumed water and cost of consumed energy is 
considered as the cost of the process. The results 
clarified that the maximum observable profit is 
obtained at the COT, Feed flow rate and steam ratio 
equal to 883 Q C, 6.02gr/min and 0.62 (gr/gr). To 
determine the effect of reactor temperature profile on 
the yield of main products and the net profit, 
optimization was carried out. The results clarified that 
energy consumption has important effect on the net 
profit. Moreover, the rate of temperature increment 
and the value of maximum temperature have 
important rolls on the maximization of olefin 
production. 
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Nomenclature 

COT coil outlet temperature (°C) 

Coef . coeficient in quadratic model 



Tfurnace 

Tskin 

XOT 

XI 

X2 

X3 
Yi 

a 

pij 

A 



temperature of reactor furnace 
(K) 

skin temperature of the reactor 
(K) 

cross over temperature (-C) 

experimental design parameter 
(feed flow rate (gr/min)) 

experimental design parameter 
(steam ratio) 

coil outlet temperature ( Q C) 

yield of products (mass %) 

constant coefficient in quadratic 
model 

coefficients in quadratic model 

experimental factor for 
predicting the relationship 
between decoking time and 
operating time (hr/hr) 
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